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In a recent publication 1 we have shown the possibility to achieve strong coupling of the quantized 
motion of a micron-sized mechanical system to the motion of a single trapped atom. In the proposed 
setup the coherent coupling between a SiN membrane and a single atom is mediated by the field 
of a high finesse cavity, and can be much larger than the relevant decoherence rates. This makes 
the well-developed tools of CQED (cavity quantum electrodynamics) with single atoms available in 
the realm of cavity optomechanics. In this paper we elaborate on this scheme and provide detailed 
derivations and technical comments. Moreover, we give numerical as well as analytical results for a 
number of possible applications for transfer of squeezed or Fock states from atom to membrane as 
well as entanglement generation, taking full account of dissipation. In the limit of strong-coupling 
the preparation and verification of non-classical states of a mesoscopic mechanical system is within 
reach. 



PACS numbers: 

I. INTRODUCTION 

The quantum regime of optomechanical systems [3J |S] 
- in particular micro or nanomechanical oscillators cou- 
pled to the optical field in a cavity - has recently received 
considerable attention, mainly owing to the experimen- 
tal progress in quantum ground state cooling [H and 
strong coupling dynamics [6j III [9l [10]. Combining 
opto-micromechanics with low-loss dielectric membranes 
pn [T^ on the one hand, with cavity QED [T3] with sin- 
gle or many atoms on the other hand, a hybrid system 
emerges that can be a testbed for experiments on coher- 
ent dynamics between microscopic (single atom or ensem- 
ble of atoms) and macroscopic (micro-mechanical oscilla- 
tor) systems. Given the already well-developed toolbox 
for the manipulation of atomic states such an interface 
can be used for indirect preparation and manipulation 



of quantum states of mesoscopic mechanical oscillators. 
Moreover, in view of applications such as quantum infor- 
mation processing, it seems timely to ask for quantum 
hybrid systems which combine the advantages of physi- 
cally different systems, each with a unique set of proper- 
ties and capabilities, in a compatible experimental setup. 
A hybrid atomic-mechanical system would be one such 
example. 

A few recent theoretical proposals advance the possi- 
bility of coupling ensembles of atoms to mechanical res- 
onators. Most generally the interaction is mediated by a 
light field that couples the mechanical resonator via the 
radiation pressure effect to either internal levels of the 
atoms [m [T21 [TB] , or to their motional degrees of freedom 
|17j . which can result e.g. in cooling of the mechanical 
resonator via a bath of atoms jH] . Also a direct coupling 
has been proposed where a magnetic tip mounted on a 
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FIG. 1: Dynamic intracavity field provides strong interface 
between the motion of a single trapped atom and the vibra- 
tions of a micron-sized membrane. 



cantilever provides a Zeeman coupling to the atomic spin 
of the Bose-condensed [T^ or ultracold atoms. A 
number of proposals discuss the possibility to couple the 
motion of a microresonator to single two level systems, 
realized e.g. in a quantum dot [21j . in a nitrogen- vacancy 
impurity in diamond |22] , or in superconducting circuits 
such as a Cooper-pair box |23l |24| , a SQUID [25l |26] or 
a flux qubit [57j. 

The direct coupling of the motion of a single micro- 
scopic body such as a single atom to a macroscopic me- 
chanical oscillator is considerably more challenging. Typ- 
ically, the interaction strength is governed by a small in- 
trinsic parameter which scales as -y/m/M ~ 10"'' — 10"'*, 
where m and M are the masses of the atom and the me- 
chanical oscillator. This is true e.g in [2E| where the 
motion of an ion in a trap is coupled to the vibrations of 
nano-electrodes providing the trap potential. An alterna- 
tive route is however possible, where an indirect cavity- 
mediated coupling circumvents the limitations imposed 
by the small mass ratio, as presented in our recent pro- 
posal Thereby a strong coupling is achievable be- 
tween a single trapped atom and the motion of a mem- 
brane, where the coupling strength can exceed the dissi- 
pative rates by a factor of ten for present or near future 
experimental parameters. 

In this article we elaborate on the mechanism described 
in our previous letter |T], and provide more details and 
applications of the scheme. The paper is structured as 
follows. Section [IT] presents an overview and qualitative 
picture of our results. In Sec. |III| the reduced master 
equation describing the cavity-mediated membrane-atom 
interaction is derived in detail, and results are presented 
in particular for the dispersive regime. Sec. |IV| specializes 
on the regime of strong membrane-atom coupling, and 
examples of state transfer are presented. In addition, 
we describe how to produce entanglement between atom 
and membrane by modulating the input laser intensity 
in time, leading to a two-mode squeezing Hamiltonian. 
Sec. |V] discusses technical details regarding the specific 
setup that we have in mind, and finally we discuss the 
result and conclude in Sec. IVII Mathematical details of 
the derivation are presented in Appendices. 



II. OVERVIEW 

In the setup proposed in pL,, the recent development 
within micromechanics with membranes in optical cav- 
ities is combined with single trapped atom cavity 
QED [13]. As shown in Fig. [T| we consider an optome- 
chanical system where a micron-sized dielectric mem- 
brane is placed in a laser driven high-finesse cavity and 
coupled through radiation pressure to the cavity field 
quadratures, with the coupling strength controlled by 
the laser power through the intracavity amplitude. In 
this setup, the membrane vibration manifests itself as 
a dynamic detuning of the driven cavity modes. For a 
cavity mode driven by a laser detuned from the cavity 
resonance this dynamic detuning translates into a dy- 
namic intracavity field intensity. If now a single atom is 
trapped in the optical dipole potential provided by the 
cavity field, the membrane vibration couples via the dy- 
namics of the optical trap to the motion of the atom, 
and vice versa. This coupling is strongly enhanced by a 
large steady-state field amplitude and the cavity finesse, 
which is a key ingredient in achieving the strong coupling 
regime. 

A. Effective Master Equation and Strong Coupling 

The focus of our analysis is a configuration where the 
cavity field serves merely as a quantum bus and can be 
effectively eliminated from the dynamics, giving rise to a 
coupled oscillator dynamics for the reduced system com- 
prising the membrane and the atom [h ~ 1], 

H = UJ^nO-i^am + ti^at alt flat " (^(flm -|- a|^) (Cat + a^t) • (1) 

In this Hamiltonian the first and second terms describe 
the bare micromechanical oscillator and harmonic motion 
of the trapped atom, respectively, with Om (flat) being the 
mechanical (atomic motion) annihilation operator, uj^ 
and Wat are the respective oscillation frequencies. The 
linear form of this interaction would provide a quantum 
interface for coherent transfer of quantum states between 
the mechanical oscillator and the atom, opening the door 
to coherent manipulation, preparation and measurement 
of micromechanical objects via well-developed tools of 
atomic physics, as will be detailed in Sec. |IV| 

However, the cavity mediated, coherent dynamics will 
compete with a number of dissipative processes, such that 
the full dynamics will be described by a master equation, 

p= -i[H,p] + {p) + Lat (p) + ic (p) • (2) 

The three Liouvillian terms describe the respective 
sources of dissipation, with L,^ including the thermal 
heating of the membrane vibration and Lat including the 
atomic momentum diffusion due to spontaneous emis- 
sion. Furthermore, a cavity-mediated coupling comes 
naturally at the price of cavity-induced decoherence via 
photon leakage, Lc- Our goal here is to construct a setup 
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FIG. 2: Linear atom-membrane coupling mediated by two 
driven cavity modes, (a) One mode is driven on the red 
side, the other on the blue side. When the mode frequen- 
cies shift due to the membrane vibration (dashed line), the 
cavity response is reduced for one mode and enhanced for 
the other, (b) Atom and membrane in equilibrium inside the 
driven cavity, (c) When the membrane vibrates around its 
equilibrium, the oppositely changing cavity response for the 
respective modes shifts the equilibrium of the combined atom 
potential. 



that the response of the cavity field amplitude to the 
atomic motion is close to maximal. Similarly, the mem- 
brane is positioned at aim half-way between a field node 
and anti-node, where the linear opto-mechanical coupling 
is maximal [TT] , see Fig. ^h) . The position Xra is chosen 
such that both fields have similar slope (with the same 
sign), thus react equally to the membrane vibration. The 
displacement of the membrane thus shifts the cavity res- 
onances, as shown by the dashed lines in Fig.[2]ja). With 
the two lasers being tuned to different sides of their re- 
spective resonances, during the membrane displacement 
one driving laser will come closer to resonance, with re- 
sulting enhanced intracavity field, and the other one far- 
ther off resonance with resulting reduced intracavity field. 
Consequently we will find that one of the atomic lattice 
potentials is getting deeper, the other one getting more 
shallow, as seen in Fig. ^c), thus shifting the atomic 
trapping potential. Due to this spatial shift being pro- 
portional to Xni, the result is an overall ~ a:;ata;m coupling 
as in Eq. ([l}. 

With this construction, the cavity field can provide 
the leverage to couple two objects with mass ratio on the 
order of 10~^^. Imagine for illustration that we were to 
achieve a similar coupling with a mechanical device like 
a seesaw: To balance the torques would require a lever 
ratio of the same order of magnitude; 15 mm on one side 
and the earth-sun distance on the other side. 



obeying the master equation (|2| with a Hamiltonian term 
([T| where the interaction between the atom and the mem- 
brane is resonant, i.e. co^ ~ Wat, and strong, i.e. the 
coupling constant G is larger than the relevant deco- 
herence rates rc,ri„,rat corresponding to the dissipative 
processes described by Lc, L^, L^t, respectively. In fact, 
we will show that for state of the art experimental pa- 



rameters small ratios (FcFi, 
reach. 



, rat)/G ~ 0.1 are within 



B. Qualitative Picture of Linear Coupling 

A strong linear coupling as described in Eq. ([T]) is 
obtained in a configuration involving two driven cavity 
modes of frequencies Wc,i and Wc,2i as shown in Fig. |2ja). 
The two modes are driven by lasers of frequencies wi and 
UJ2, respectively, where the first laser is tuned to the red 
side of its cavity resonance, uji —ujd < 0, and the second 
laser is tuned to the blue side, 1^2 ~<^c,2 > 0- By a proper 
choice of cavity modes, and with an internal structure of 
the specific atom as shown in Fig. |3] both lasers sepa- 
rately provide red-detuned optical lattices, which com- 
bine into a potential where we trap a single atom in one 
of the wells, see Fig. |2];b). With wave vectors fci ^ ^2 
and assuming equally large intracavity amplitudes, the 
two lattices have opposite slopes at the equilibrium po- 
sition Xat of the atom. Moreover, the particular well is 
chosen such that the slopes are close to maximal, such 



III. MODEL FOR CAVITY MEDIATED 
MEMBRANE-ATOM COUPLING 

After the qualitative description in the last section, we 
move on to a detailed presentation of the system consist- 
ing of a moving atom and a vibrating membrane coupled 
to driven cavity modes. Further, we will show how to 
obtain the reduced atom-membrane dynamics described 
by Eq. (2^ by eliminating the cavity degrees of freedom, 
and finally we will identify the regime of strong coupling. 



A. Detailed Derivation of Effective Master 
Equation 

1. Full master equation 

Our starting point is the complete master equation for 
the density operator W describing the dynamics of cavity 
modes, atom and membrane motion. 



{W), (3) 



with the coherent dynamics contained in the system 
Hamiltonian H^^^, 



Hsys — -^motion ~t~ Hq 
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FIG. 3: Two lasers with frequencies wi and UI2 respectively 
drive two different internal atomic transitions with detunings 
5. 



Here i?niotion takes into account the free harmonic motion 
of the membrane, modeled as a single-mode oscillator, 
and the kinetic energy of the atom with momentum Pat , 



(4) 



Further, contains the free cavity Hamiltonian, as well 
as the effect of the atomic motion and the membrane 
vibration on the cavity field. We will postpone its dis- 
cussion to the next section, where we give the concrete 
form of He for various setups, and here first address the 
remaining terms in the master equation ([s]). 

The Lindblad terms in the master equation ([s]) describe 
dissipation of the membrane (Lm), the atom (Lat) and 
the cavity modes {Li) respectively. Here i labels the cav- 
ity modes with photon annihilation operator Ai obeying 
= 5ij. Cavity decay at an amplitude decay rate 
Kj is described by, 

U{W) = K,V[A,]{W), 

where we use the shorthand notation 

V[a]{W) = 2aWa^ ~ {a} a, W}+ 

for a Lindblad term with jump operator a. 

Lat describes the noise processes acting purely on the 
atoms, and could also include controlled dissipation such 
as Raman cooling. Here we assume each driven mode 
with wave number ki to couple to a different atomic tran- 
sition |0) ^ as sketched in Fig. [s] and focus on the 
photon recoil during spontaneous emission from the ex- 
cited state |i) with rate 7^ to the common ground-state 
|0). The effect of the photon recoil on the atomic motion 
is described by the following Lindblad term. 



2 / d7?S,(i?)e^'"='^-\A^(^M^e-''''='^-y^i;(^ 



{M,(Xat),W^}_ 



(5) 



where 7^ is the spontaneous emission rate and Si is the 
saturation parameter for transition i [331. Si(fi) are 



(6) 



(even) geometric functions whose exact expressions de- 
pend on the chosen transitions, and Ui(x^\) describe the 
spatial intensity profiles of the cavity modes with x^t the 
atomic position in the cavity. In the next section we 
will discuss this term in more detail in the Lamb-Dicke 
regime, where it simplifies considerably. 

Finally, Lm describes the membrane thermal contact 
via the finite temperature suspension, modeled as inter- 
action with a thermal bath, 

with 7ni the natural linewidth of the mechanical reso- 
nance, and fini its mean occupation in thermal equilib- 
rium. The heating rate Fm = 7m^m — ksT /{TlQ^) is 
related to the temperature T of the contact and the me- 
chanical quality factor In addition to the thermal 
contact, we include membrane heating due to absorption 
of laser power. In fact, a fairly cautious estimate detailed 
in section |V] shows that with standard cryogenic precool- 
ing the natural lower limit for the temperature T is set 
by light absorption within the membrane. 

The model presented so far makes the following as- 
sumptions: (1) atomic motion is accurately described by 
a ID model, with the transverse confinement provided by 
the Gaussian intensity profile of the cavity fields, (2) in- 
ternal atomic dynamics can be eliminated, assuming the 
laser drive to be sufficiently detuned from the atomic res- 
onances (cf. Fig. |3|, and (3) negligible internal coupling 
of the chosen membrane mode to vibrations of higher 
energy, allowing a single mode approximation. 

2. Linearization Around Equilibrium 

We will now proceed to discuss the cavity and inter- 
action Hamiltonian He- Note first that despite the very 
different physical nature of the atom and the membrane, 
their effect on the cavity field can be collected in a unified 
description, where the cavity modes Ai with frequencies 
t^ci see an index of refraction which depends on the re- 
spective positions Xat and of the atom and the mem- 
brane along the cavity axis. This description assumes a 
Born-Oppenheimer type approximation, where the slow 
atom/membrane motion compared to the optical cavity 
frequencies allows a separation of timescales. 

Single Mode Setup: We now first want to consider a 
setup with only a single driven cavity mode, i = 1, in 
order to illustrate a number of conceptual points. The 
generalization to the two-mode setup is then immediate. 
For a single mode the cavity Hamiltonian in the mas- 
ter equation (|3|, taken in a frame rotating with the laser 
frequency u; (dropping the index i), is 

H, - [u, (xat, x„,) - Lo] A^A + E (e^-^At + h.c.) . (7) 

The first term is the cavity free energy in the rotating 
frame, which depends parametrically on the atom (mem- 
brane) position. The second term describes the laser 
drive of power P, such that E = yj2PK/fiu}c. 
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The strong drive field creates a steady-state intracav- 
ity field with amplitude a ^ 1, which in turn provides 
a trap potential for the atom at a certain equilibrium 
point Xat and mean force on the membrane, displacing 
it to a slightly shifted position Xm- We are interested 
in the dynamics of the fluctuations of cavity amplitude 
and atom/membrane position around these equilibrium 
values. It is therefore convenient to move to a displaced 
frame where the dynamics is described by the fluctua- 
tions a around the steady-state field, 



The term in the second line of Eq. (12) describes 



A = 



a, 



(8) 



and the fiuctuations Sxi^t said Sx^ around the equilibrium 
atom and membrane positions. 



Sx„ 



(9) 



Along this line, we expand the cavity mode frequency 
LUc {xnt, Xni) around steady-state. 



UJc {Xat,Xni) 

1 

+ 2 L'^at 



1 



(10) 



[dl 



+ [l9at,m'^c] 5x.^t5x 



with dgx and short for the partial derivative with re- 
spect to a; at and Xm, evaluated at the atom and mem- 
brane equilibrium points. When the expansions (|8]),([9| 
and ( 10 1 are used in the full master equation ([s]) with 
Hamiltonian ([t]) , the steady state amplitude a and equi- 
librium positions of atom and membrane can be deter- 
mined self-consistently by demanding that all terms van- 
ish which are linear in fluctuation operators 5xat^ 5xa\ 
and a. In particular one finds for the intra-cavity ampli- 
tude (see [22 for details). 



(11) 



{lO — (jjj?) -I- IK 

The laser phase can be chosen for convenience such as 
to make a real. 

In the resulting Hamiltonian all linear terms are thus 
systematically removed and the dynamics is governed by 
an effective Hamiltonian 



at.m 



• {ljI — Ljj)a} a 



5x1 



+ a [9mWc] <^a:^in(a + a^) + ^ [d^^^c 
+ f [^at^^c] 5xl,{a + a^), 



[dl 



5xl^{a- 



5x1^ 

f a^) 
(12) 



where terms of fourth order in fluctuations (zeroth order 
in cavity amplitude) have been neglected. 

In the first line we find the free energy of the cavity and 
the optical potential for the atom, providing a harmonic 
trap with a frequency determined by mw^^ = [cJ^^Wc] . 
The corresponding term for the membrane provides a 
small correction to its mechanical frequency and can be 
neglected. 



a direct linear atom-membrane coupling of the form 

~ 5diroct(aat +alt)(«m + a3n)i where feat = 4t(aat + alt) 
and Sxya — ^m(im + '2m)- The zero point fiuctuations are 
given by. 



2ma;at 



1/2 



\2MuJn 



1/2 



and now refers to the shifted frame for the membrane. 
Assuming that the cavity field provides the atomic trap 
as discussed above, with a trap frequency close to that 
of the membrane vibration, Wat ~ ^m, it can be checked 
easily that the direct coupling will be hampered by the 
small mass ratio. 



^?dircct 



(^m/4t)Wat ~ yJm/Mu}^ 



thus would not reach the strong coupling regime for a 
single atom. We will see that the cavity-mediated, in- 
direct coupling can be many orders of magnitude larger, 
such that the direct coupling can be safely neglected in 
the following. 



The third line in Eq. ( 12 ) describes a membrane-cavity 
interaction. As was discussed in detail in [30] a proper 
choice of the membrane position along the cavity axis can 
make either the first or the second term dominant. In the 
later case the cavity field couples to Sx^ ~ (am + aJn)^; 
which has interesting applications for measuring occupa- 
tion numbers of the membrane. However this term is 
typically rather small as it scales like d^^LOct^ ~ (^c^m)^ 
and is thus of second order in the corresponding Lamb- 
Dicke parameter. In the following we will neglect this 
second order term and keep only the first one, where the 
cavity couples linearly to the membrane fiuctuations. 

The linear-coupling term for the atom vanishes, as for 
a single cavity mode the atomic equilibrium position is 
defined by [Satij-'c] — 0. Thus, only the quadratic, para- 
metric term given in the last line of Eq. ( 12 1 contributes 
to the atom-cavity coupling. It is perfectly possible to 
proceed from here and to derive an effective coupling of 
the atom to the membrane. However, this coupling will 
be ~ Xa^x'^^ and thus not of the desired form given in 
Eq. 0. 

Two Mode Setup: Creating a linear atom-cavity cou- 
pling requires non-vanishing cavity field slopes at the 
mean position of the atom, [(Jaf^c.i] 7^ 0. To this end 
one has to require an external trap for the atom shifting 
it away from a lattice extremum. An elegant alternative 
is to use two driven cavity modes {i — 1,2) with the 
atomic equilibrium position at an extremum of the com- 
bined optical potential, and at the same time at a point 
of maximal slope of the individual cavity fields. Let us 
therefore study the cavity Hamiltonian in detail for 
this case of two driven cavity modes. In a frame rotating 
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with the laser frequencies uji we have 

1=1.2 

+ ^£;,(e*^-4 + h.c.), (13) 

with mode frequencies Wc.i [xaXiXm] given by, 

Wc,i (a;at,a;m) = " [5o,i/^m]a;m + U^u^ix-^t)- (14) 



The second term in ( 14 1 describes the dynamic cavity 



detuning due to vibrational fluctuations of the thin di- 
electric membrane, with single-photon coupling jlll I30j 



go,. = (CAX, 



(15) 



Here L is the cavity length and — 2r sin(2A:ia;m)/[l — 
cos^(2fciXni)]^^^ is a correction factor which takes into 
account the finite amplitude reflectivity r of the mem- 
brane, as well as the distance to the cavity field node 
where the field is zero and thus insensitive to the mem- 
brane motion. Note that the special case = 1 is famil- 
iar from optomechanics with a perfectly refiecting moving 
mirror. By a proper choice of membrane location Xm it 
is possible to achieve ~ 2r for both fields. 



The third term in ( 14 ) describes how the driven optical 



modes provide a lattice potential for the atom along the 
cavity axis, with the spatial intensity profile 

Ui(x) — sin^(fcix), 

and a lattice potential strength determined by the AC 
Stark shift (per photon). 



(16) 



Here fio is the vacuum Rabi frequency and b is the detun- 
ing of the lasers from the respective atomic transitions, 
assumed equal for simplicity (see Fig. |3|. 

The linearization of the dynamics around the equilib- 
rium mean values is done as for the single-mode case 
discussed previously. The intracavity amphtudes are 
ai = Eie^'^'/{Ai + iKi) with Ei the drive strength of 
mode i, and the phase is chosen to make real. 



A,; = 



is the laser detuning relative to the cav- 



ity mode. The following derivation is in principle gen- 
eral regarding the number of driven modes, the mode 
parameters ai,Ki and the optomechanical coupling go.i- 
Without loss of generality we will assume in the follow- 
ing a symmetric two-mode case with ai = a, Ki = k and 
9o,z = 90- 

The expansion is again very similar to the setup for 
a single mode. The main differences concern the atomic 
degrees of freedom. The atomic mean position x^t is 
determined by vanishing first derivative of the total field. 



and the trap frequency Wat of the harmonic potential at 
this position is accordingly 

mtol = Uoa^klCix^t), Ci^) = u"{x)/kl 

Most notably the ac Stark shift term Uq J2i'Ui{x at) A^ A 
gives rise also to a linear atom-cavity interaction, as the 
individual terms ~ u'j^^x^t) can be nonzero despite the 
condition on vanishing derivative of the total field ([Tt]), 



Ul(2) (a^at) ^ Ui(2){Xat) ± V^{Xat){aat + 4t) (18) 

with 0(x) — u'i{x)/ki and Lamb-Dicke parameter 77 = 
fci^at- A significant slope 9{xat) and hence significant 
coupling is achieved for two modes by a careful choice of 
atomic site within the cavity, far from extremum points 
of the individual lattice modes, as we discuss in Sec. |Vj 
In this way it is straightforward to expand and lin- 
earize the cavity Hamiltonian for the two mode setup in 
Eq. ( 13 1. When combined with the kinetic energy of the 



atom this results overall in a linearized Hamiltonian 77=- 



of the full master equation in Eq. ([3| 



int 



(19) 



with a free energy 



and linear membrane-cavity and atom-cavity interaction 

-ffint = .gm(am + ajn)[(ai + «l) + («2 + 4)] 

+ gat(aat + alt)[(ai + 4) ~ (^2 + 4)] 
with coupling strengths. 



ffat = Uoarie{xat), 



9n 



9oa- 



(20) 



In Hgys all parametric coupling terms have been ne- 
glected, as they will be smaller by the atomic Lamb-Dicke 
factor r/ or by the much smaller Lamb-Dicke factor cor- 
responding to the membrane motion, as discussed previ- 
ously. In the limit of large cavity amplitude we also drop 
all terms of zeroth order in a. For later use it will be 
convenient to reexpress the interaction in the form 



[' 



(a* + 4). 



with operators Fi describing the forces exerted by the 
atom and membrane motion on the cavity fields, 

9=4^K+9if (21) 



Fi.2= -—a, 
9 



I 5at 
± flat 



"'(a^at) = 0, u{x) ^ Ui{x) + U2{x) . (17) of 



Before we derive the cavity-mediated atom-membrane 
coupling, we will finally discuss the atomic Lindblad term 
Eq. (|5]) . In a Lamb-Dicke expansion around the atomic 
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equilibrium position in the optical potential this Lind- 
blad term takes the form of a momentum diffusion master 
equation 

Ut{W) = ^-fV[a^, + al]{W), 
with a diffusion rate, 

Tat = ri'sel [2 - (4/5)u(a;at)] (22) 

where the saturation parameter i s no w explicitly given by 
Se = [a^lo/S]'^. The expression (22 1 in the end depends 



on the particular atomic transition and the specific ge- 
ometry; the factor (4/5) is specific for transitions with 
Am = but for other transitions it is still of order unity. 

Let us remark that in fact it is possible to solve the 
master equation of the full system exactly (e.g. by means 
of the methods given in Appendix |b]) , and indeed there 
can be rich physics to be explored in the regimes not 
considered here. However, the focus of this paper is 
the regime where the cavity modes can be eliminated, 
5at,5m ^ max{K, A}, which is not only more relevant 
from an experimental point of view, but also allows for 
analytical, transparent results which highlight the phys- 
ical properties of the system. 



3. Adiabatic Elimination of Cavity Field and Effective 
Master Equation 

We are now in the position to derive an effective cou- 
pling mediated by the cavity modes. The idea is to use 
a parameter regime where the cavity dynamics is essen- 
tially unperturbed by the motion of the membrane and 
the atom, and solely mediates interaction between the 
two. The corresponding requirement is fast cavity dy- 
namics, g <ti K or g ^ \Ai ± lu^\. For optomechan- 
ical cooling the former condition is the more common 
requirement, but since the resulting strong dissipation 
through the cavity decay would harm the coherent cavity- 
mediated dynamics, we choose a regime where A^ are the 
large parameters. Here fluctuations in the cavity quadra- 
tures are fast variables which adiabatically follow the dy- 
namics of the position fluctuations of the atom and the 
membrane. In order to achieve strong interaction we fur- 
ther assume atom and membrane to be on resonance. 



Wat 



(23) 



The formal procedure for eliminating the optical modes, 
as described in detail in Appendix [X] is to perform adia- 
batic elimination using standard techniques 133] . We find 



that the linearized atom-membrane-cavity dynamics ( 19 ) 
gives rise to the effective master equation ^ with 



H = Winamam + Watflat^at + ^^at-m 



The last term iJat-^m (A3 1 represents the cavity-induced 
atom-membrane coupling and a correction to the free 



motion, and can be extracted from the coherent part of 



the cavity-mediated Liouvillian Lc-mcd ( Al ). In detail it 
reads. 



Fa 



E 



K + i{Ai - Win) 



F, 



K + i{Ai + LUn-i) 



F 



h.c. 



(24) 



The cavity decay translates into correlated decay Lc{p) 



( A4| for atom and membrane, where in the rotating wave 
approximation (RWA) each optical mode i contributes 
cooling (I?[Ji]) and heating (I?[f/]) associated with emis- 
sion of sideband photons at either side of the driving 
laser, that is, at one of the two frequencies LOi ± Wm, 



g^K 



k2 -I- (A, - Wm)^ 



T^iFhip) 



(25) 



An emission event is accompanied by the creation or an- 
nihilation of a quantum in either the atomic motion or 
the membrane vibration. For a near resonant system 
(wm — Wat) these two possibilities are indistinguishable, 
such that both processes happen in a coherent fashion. 
Therefore, the jump operators Fi are linear combinations 
of the corresponding annihilation operators Oat and a^. 



4-. Coupling in the Dispersive Regime 

So far, we have derived expressions for the cavity- 
mediated interaction (24 1 as well as its inevitable com- 
panion, dissipation through the cavity decay (25 1. The 



remaining challenge is to reach the strong coupling 
regime for the reduced system, with effective coupling 
strength G which is much larger than all decay rates, 
G ^ Tc, Fin, Fat- Let us first consider the relation to the 
cavity-induced dissipation described by Lc- 

Our first observation is that the atom-membrane cou- 



pling (24 1 is maximized for equal and opposite detunings. 



for which the two cavity modes respond equally and op- 
positely to the membrane vibration. Evaluating the ef- 
fective Hamiltonian (pil for this special case we find 



= -G 



(a,n + a|n)(aat + alt) 



-I- ie (^amflat - a]jja^t 



where we dropped a global energy shift. The effective 
coupling strength G is given by 



G^ 



2gni5at(A - t^m) , 2g,ngat ( A -|- W,n) 



(A - 



(A + CJn 



8 



From the observation that the rate of cavity induced 
decoherence in Lc, see Eq. (25), scales like ^ 1/A^ 
whereas the cavity mediated interaction G ~ 1/^; 
draw the conclusion that the dispersive limit is natu- 
ral for suppressing dissipation. Focusing on the regime 
where |A| is the largest parameter, |A| ^ Wm, n, the cor- 



rection e to a pure (a„ 
neghgible. 



7t 



)(^ 



^)-interaction is 



A2 + ^2 _ ^2 



< 1. 



Thus the coherent dynamics in the reduced master equa- 
tion ^ is effectively given by the Hamiltonian _ff in ([T|). 
This is the main result of our investigation. To zeroth or- 
der in k/A, ojm/A the coupling constant G has the simple 
form, 



G 



45m5at 

A ' 



Regarding the cavity-induced decoherence processes de- 
scribed by Eq. (25 1, the combination of a red-detuned 
(Ai = A < 0) and a blue-detuned (A2 = — A) laser drive 
can be interpreted as simultaneous cooling and heating 
processes. The rate of cooling r+ via mode 1 equals the 
rate of heating via mode 2, and vice versa with rate F", 



proposed setup. In particular we discuss how to obtain 
a time- dependent atom-membrane coupling G{t), which 
can be advantageous e.g. for entanglement creation as 
discussed in section HVl 



1. Single driven mode combined with external trap 

As briefly mentioned previously, as an alternative to 
the two-mode setup one could use a single driven mode 
combined with an external atom trap. This trap would 
shift the atom away from the lattice extremum, to an 
equilibrium point Xat where the cavity field has finite 
slope, u[{xi^t) 7^ 0, which can be significant for a prop- 
erly chosen atom location within the cavity. With cou- 
pling only to a single mode, the force Fi is given by 
(21 1 whereas F2 — 0. Consequently, with only one cav- 



ity mode mediating the membrane vibration, the largest 
displacement of the atomic mean position is only half as 
large compared to the case when two cavity modes are 
shifted out of phase. Hence the resulting coupling con- 
stant is only half of its maximum value. 



single— mode 



2ffmgat 

A ■ 



L,{p)=^(v[F,]ip)+V[F^]ip) 



V[F}]{p) + V[F,]{p)), (26) 



with the rates given by. 



F^" 



(ffm + git) 

(A±W„i)2 



In our attempt to minimize dissipation we additionally 
note that the relation between the coupling constants 
(7at and 5ni is of importance. The ratio of dissipation to 
coupling strength is proportional to, 



F±/Gcx 



gm + ffat 
.9mgat 



This implies that the mediated atom-membrane inter- 
action is most efficient when the two oscillators couple 
equally strongly to the cavity modes, = g^t- Un- 
der this condition and to lowest order in K/AjWm/A the 
cooling / heating rates are in fact equal. 



F^ 



Fe-G 



A' 



(27) 



a factor k/A ^ 1 smaller than the coupling constant G. 



B. Alternative setups 

In this section we extend the previous discussion to 
jive a hint about alternative mode configurations for the 



2. Time- dependent coupling constant G{t) 

In principle, it is possible to achieve a time-dependent 
coupling constant G{t) by modulating the laser power, re- 
sulting in a time-dependent intracavity amplitude a{t). 
However, the basic problem with simply introducing 
time-dependent a{t) in the two-mode setup, is that the 
atom potential will be time-dependent as well, which 
could heat up the atomic motion. Therefore we consider 
the following modified setup, where we either use an ex- 
ternal trap for the atom as above, or use a second mode 
mainly to provide the trapping potential. Either way, the 
role of the first mode is to mediate the atom-membrane 
interaction with modulated strength, ai{t) = aic{t). 

Without going into details, the idea of the two-mode 
case is to drive the second mode such that the corre- 
sponding intra-cavity field becomes very strong, a2 ^ 
ai, with Qfi the amplitude of the coupling mode. Due to 
the very small ratio a\l ^ 1, the first mode hardly in- 
fluences the atomic potential at all, and the mean atom 
position Xat is given by vl^ixzx) = 0. Furthermore the 
atomic frequency is determined by the curvature of the 
second fleld, mt^at = Q^2'"2 (^^at)- Since M2 = at the 
atomic equilibrium point, the second mode will not con- 
tribute to the linear atom-cavity coupling. 

With the atom-membrane interaction mediated by 
only a single mode, as discussed above, the couphng con- 
stant will be half of the maximum value. We find the 
resulting time-dependent atom-membrane coupling. 



G 



two— mode 



2gmffat 2 



A 



c\t). 
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We will come back to this possibility of making the cou- 
pling explicitly time dependent in our discussion of co- 
herent evolution, in particular of a protocol to gener- 
ate entangled states of the atom and the membrane, see 
Sec. HVEl 

We will remark on yet another type of setup in the 
outlook of this paper (Sec. VI), namely how to imple- 



ment an optomechanical Jaynes-Cummings model where 
atomic internal degrees of freedom are used instead of 
the atomic motion. Together with the examples of this 
section, the discussion illustrates that the present setup 
actually provides a toolbox for engineering various inter- 
actions and different types of dynamics. 



IV. COHERENT EVOLUTION IN THE 
STRONG COUPLING REGIME 

In the previous section it was shown that we can im- 
plement a linear atom-membrane interaction (jlj with 
the proposed two-mode setup operated in the dispersive 
regime, and that this interaction can be fast on the time 
scale of relevant decoherence rates in this system. In this 
section we will study a few applications, which become 
accessible in this regime. 

Note first that the coherent evolution governed by this 
Hamiltonian transfers a state from the atom to the mem- 
brane, and vice versa, in a time given by 

ts - V(2G), 

such that |V'i)at|i/'2)m |'02)at|?Ai)m, up to local ro- 
tations. The state swap mechanism appears naturally 
in the interaction picture; for resonant coupling = 
Wm ^ G the Hamiltonian takes a beam-splitter form in 
the RWA, 

Hi~G (flataJn + h.c.) . 

Particularly intriguing is the ability to use the state 
transfer to control the mechanical state through the avail- 
able atomic physics toolbox. However, as already dis- 
cussed, the coherent interaction is accompanied by sev- 
eral sources of noise which in the end reduce the fidelity 
of the state transfer. Strong coupling is therefore estab- 
lished by fulfilling, additionally to the resonance condi- 
tion, the following set of conditions, 

G ^ Fat, Fm, Fc- 



(28) 

In section IV Al we will summarize and comment on the 
optimization of parameters which is necessary in order 
to reach the strong coupling regime, following [1]. Here 
we illustrate the strong coupling in the presence of noise 
with three specific examples of state transfer from atom 
to membrane: coherent and squeezed state as well as a 
Fock state. Aiming at a clear picture of the effect of 
dissipation, we assume all dissipation rates to be equally 
strong and define a ratio /, 



Note that we here put the membrane heating and cool- 
ing rates equal, (nm -I- l)7m — n-mlm = Fm, and as- 
sume the dispersive regime where F^ — > F^. Our aim 
is to find the acceptable noise level which still allows for 
state transfer, by solving the master equation ^ with 
the Hamiltonian ([T|) exactly. These numerical solutions 
are combined with analytical calculations based on the 
RWA as described above. In fact, for Gaussian states it 
is straightforward to analytically solve for the time evo- 
lution; the details of the derivation are presented in Ap- 
pendix [B] Interestingly, with the generalized technique 
presented in Appendix |B 2[ the impact of noise on the 
evolution of Non-Gaussian states, e.g. Fock states can 
be derived as well. In all three examples the noise intro- 
duces a thermal population rig during the time-interval 
is needed for state transfer. We find. 



Another interesting application is cooling of the mem- 
brane via coherent state swap, as will be presented in 
subsection |IVD[ Finally in IV E we present a way to en- 
tangle atom and membrane, using a time-dependent cou- 
pling G{t) - as discussed in Sec. IIIB 2 - which enhances 



exactly those terms which are neglected in the RWA. 



A. Coherent state swap 

Our first example is the transfer of a coherent state 
1/3) from the atom to the membrane. The perfect state 
swap evolves a state |0)m|/3)at into |/3e*'^^*^)m|0)at with 
the phase <f>(t) governed by the system Hamiltonian. Here 
we use the fidelity F, defined as the overlap between the 
original atomic wave function and the final membrane 
wave function, as a figure of merit for the effect of noise 
during the state transfer. Figures [4]^a,b) show the fidelity 
for transfer of two different coherent states, \(3 = 1) and 
1/3 = 5). Particularly interesting is the fidelity of the state 
swap, i.e. for t = tg where F is close to maximal but still 
deteriorated due to dissipation. In fact, the dependence 
of the state swap fidelity (at t = ts) on the noise ratio / 
follows the analytical result ( |B4[ ) derived in Appendix [b] 
in the RWA, 



1 



l + 7:f 



(30) 



This simple analytical result, which very well matches 
the exact numerical solution shown in Fig. |4|c), states 
that with noise levels below 10% we can expect a state 
transfer fidelity above 75%. 



B. Squeezed state transfer 



/ 



Fc 
G 



Fat 

~G 



I\n 

G 



The second example is the transfer of an atomic 
(29) squeezed state |^) with minimal variance AX^^ = (l/2)s; 
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FIG. 4; Fidelity for transfer of coherent state |/3) from atom 
to membrane. (a,b) Fidelity as function of time for transfer 
of a state with (a) /3 = 1 and (b) /3 = 5 for various values 
of the dissipation ratio /, with fixed G/cOm = 0.034. Here 
/ = 0.01 (black solid line), / = 0.05 (orange dotted line) 
and / = 0.10 (blue dashed line). The little wiggles are due 
to counter-rotating terms, (c) Snapshot at t = Tr/2G: fidelity 
for transfer of state with /3 = 1 as a function of the dissipation 
ratio /. 



here s < 1 denotes a state squeezed along the X quadra- 
ture. Such a state can be constructed using for ex- 
ample the parametric coupling to the cavity field, 



,t ^2 



(fli -|- a|) which was briefly mentioned in sec- 



tion |IIT] Ideally the swap operation transfers the atomic 
minimal variance to the membrane state, |/3)in|C)at — *■ 
|^e'''')m|/?e"^)at- Dissipation however broadens the vari- 
ance during the swap operation. 



with s{ts) > s(0). Fig. [5]^ a) shows the minimal variance 
of the membrane, reaching its lowest value after half a 
period {t = tg) when the squeezed state has been trans- 
ferred from the atom to the membrane. Obviously larger 
dissipation ratio / results in less squeezing transferred. 
In Appendix [B] we derive the following analytical expres- 
sion in the RWA for the dependence of the squeezing 
parameter s on the dissipation rate /, 




t=7T/G 
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FIG. 5: Transfer of squeezed atom state with initial variance 
AXft = /2 to the membrane, (a) Snapshot of Wigner 
functions (upper row - atom, lower row - membrane) for 
/ = 0.05 at t = 0, i = 7r/2G and t = tv/G. (b) Minimum 
membrane variance as function of time for different dissipa- 
tion ratios f ~ (black solid line), / — 0.05 (orange dot- 
ted line) and / = 0.10 (blue dashed line), (c) Transferred 
squeezing (in dB; S{t) — 10 log\Q[s(t)](dB)) as function of 
the dissipation ratio /, for fixed G/um ~ 0.034. The initial 
atom squeezing is given by s(0) = (purple dashed line) 
corresponding to —8.7 dB, and s(0) — e^^ (green solid line) 
corresponding to —4.3 dB, respectively. 



s iU) = s(0) + 2nf. 

Fig. jsjb) shows snapshots of the atom and membrane 
Wigner functions dXt — Q,t — ts and t = 2ts] one clearly 
sees how the dissipation broadens the variances. Fur- 
thermore, due to the coherent evolution, the squeezed 
membrane quadrature X'^ is not necessarily equal to the 
squeezed atom quadrature Xat. 

In Fig. [sjjc) we show how the membrane minimal vari- 
ance increases with the noise ratio / for two specific ex- 
amples of initial minimal variance of the atom. The exact 
result confirms the loss of squeezing given by the expres- 
sion above for s{ts). With noise levels below / ~ 10% an 
initial atom squeezing below —4.3 dB allows for squeezing 
of the membrane. 



C. Fock state transfer 

The previous two sections dealt with the engineering 
of Gaussian states of the mechanical resonator, while the 
ultimate goal would of course be to apply these meth- 
ods to create more non-classical state, e.g. states with 
negative Wigner functions. 

As a last example we therefore present the transfer of 
a Fock state with n = 1 from the atom to the membrane. 
Assuming the membrane to be ground-state cooled, the 
ideal evolution reads |0)m|'T' — l)at |n = l)m|0)at- 
The quantum properties of the Fock state are best illus- 
trated by the negative value of its Wigner function at the 
origin. Fig. [6] shows cuts through the Wigner function 
for three instants in time, t = 0, t = tg and t = 2ts. 
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Thermalization during the state transfer is reflected in 
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FIG. 6: Transfer of the Fock state |7i = 1) from atom to 
ground-state prepared membrane. Snapshot cuts through the 
Wigner functions at times t — 0, t = n/2G and t — n/G, for 
fixed dissipation ratio / = 0.05 and G/uJm = 0.034. 

the decreasing Wigner function negativity for each state 
swap. A convenient flgure of merit for the thermalization 
is therefore the value of the membrane Wigner function 
Wmifi, f3* ,t) at /3 = /3* = relative to the corresponding 
(absolute) value for a Fock state u'f(0,0), 



^w(t) 



,(0,0, t) 



kF(0,0)| 



An analytic expression for N^{t) in the RWA is derived 
in Appendix [b| see Eq. (BlOl. In Fig. [t] we present a 
case of particular interest, namely the membrane Wigner 
function negativity after the first state swap (t — tg), 
which depends on the dissipation ratio / according to. 



-l + 27r/ 



The quantum properties of the Fock state \n — 1) are 
transferred to the membrane if the dissipation ratio is 
sufficiently small, 27r/ < 1. For an experimentally fea- 
sible noise ratio / = 0.1, 14% of the Wigner function 
negativity is preserved during the swap operation. 



D. Membrane cooling through state swap 

With the present setup, we see two routes towards 
preparing the membrane ground state. The first route 
is along the lines of cavity cooling; for example cool- 
ing the membrane via the cavity decay, or via an ex- 
ternally controlled Raman atom cooling with rate Fr in 
the ground state cooling regime G <C Fr, <C cOm- The sec- 
ond route is to perform a state swap and hence transfer 
the ground state to the membrane from the (previously 




FIG. 7: Analytical result for the relative membrane Wigner 
function negativity A'^w at time t = ta as a, function of the 
dissipation ratio /. 



cooled) atom. Comparing the two routes, we find that 
the effective rate for state swap is much higher than for 
coohng, Fswap G » Tcooi ^ G^/Fr , and that the state 
swap leads to a final occupation which is a factor G/Fr 
lower than for cooling. 



-swap 



/ G 



-cool 



G \2uj„ 



Note that the expression for n™^^^ can be optimized with 
respect to G, since a large coupling strength on one hand 
decreases the noise level / but on the other hand increases 
the residual final occupation in the second term. Replac- 
ing / with Fc/G and considering as a fixed parameter, 
one obtains the minimum fij^'^P for Gopt = \/2ttT\mj^. 

Concluding that state swap cooling is more efficient 
than indirect Raman cooling, it is still interesting to make 
a comparison with typical cavity cooling in the good cav- 
ity regime where now k <§; uj^ and in the perturbative 
regime where f/m ^ k. In this case, the cooling rate 
scales as Fc ~ 5m/'* the final occupancy is 



Fiji't 

„2 



K 

2w„ 



We conclude that nf"^^ < fif for large enough atom- 
membrane coupling, G/tt > g^l K. 

One drawback with the state swap procedure is that it 
only works for a precooled membrane; the anharmonic- 
ity of the atom well supports the transfer of only a few 
quanta from the membrane, say riwcii ^5—10. The sit- 
uation looks better in a generalized setup with A^ atoms 
distributed over the lattice. In the ideal case of no atom- 
atom interaction and identical atom site conditions, the 
effective coupling G is enhanced by a factor y/N and 
the single atom operator Cat can be substituted by the 
center-of-mass operator — {^/VN) J2^=i ^atj- Here 
we can consider introducing even further anharmonicity 
of the atomic wells to prevent an individual atom to be 
multiply excited, in which case the center-of-mass mode 
can support the transfer and storage of N excitations 
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from the hot membrane, thus allowing the swap of a fairly 
large thermal membrane occupation. 



E. Entanglement 

In this subsection we lay out the prospects of observing 
entanglement between the atom and the membrane. The 
major obstacle in this regard is the coupling of the system 
to the environment. Even when we assume the mem- 
brane to be prepared in the ground-state initially, the 
system will quickly heat up and entanglement is lost, at 
least if one just considers the usual static coupling. Here 
we point out, instead, a method for generating entangle- 
ment based on a time-dependent modulation of the input 
laser intensity that controls the atom-membrane interac- 
tion strength. In fact, this scheme can be employed gen- 
erally in optomechanically coupled mechanical systems. 
The scheme turns out to be relatively robust against the 
impact of the dissipation channels. By modulating the 
atom-membrane coupling strength in time one can re- 
alize a non-degenerate parametric amplifier (two-mode 
squeezing) which induces strong quantum correlations 
between atom and membrane despite the simultaneously 
occuring heating of the system. We consider the linear 
membrane-atom interaction Eq. ([T]) with the coupling 
constant modulated according to. 



'Wat 



Git) = Gcos'^{Ldt), Ld^- 

In order to allow for a modulation of the coupling 
strength without modulating the trapping frequency the 
setup needs to be modified as discussed in Sec. |III| B. 
Switching into the interaction picture, we find that in 
contrast to the case of constant coupling G previously 
discussed, the coupling term here effectively transforms 
into a parametric amplifier part and a slowly oscillating 
beam splitter part (in RWA) , 

Hi ~ ^[a^aat + h.c] + ^[at„a,te'(-"— )* + h.c], 

and contributions that are oscillating fast with respect to 
the time-scale of G, and hence have negligible influence. 
It is the parametric amplifier part that can be exploited 
to generate strong correlations. 

As a measure of entanglement we employ the logarith- 
mic negativity [Mj [351 [SB] . For a Gaussian state it can 
be computed directly from the elements of the covari- 
ance matrix [36^, whose time evolution was derived in 
Appendix B. 

Fig. ([s]) displays the generation of entanglement (a) 
and the increase of the atomic excitation number riat — 
(a^tOat) (b). This number should not exceed a threshold 
value of Tiat ~ 5 — 10 in order to keep the effects of the 
anharmonicity of the trap negligible. Note that we chose 
a relatively large difference in the oscillation frequencies, 
"^at/wni = 1.1, in order to suppress the infiuence of the 
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FIG. 8: Generation of entanglement between atom and mem- 
brane by modulating the coupling strength, (a) Logarithmic 
negativity as a function of time for diflerent values of the 
dissipation ratio / = 0.1 (blue curve), / — 0.05 (red curve) 
and / — 0.01 (black curve), (b) Corresponding increase of 
the atomic excitation number Wat ~ (a^taat). For these plots 
we chose G/wm = 0.034 and i^at/ti'm = 1.1 and assume the 
membrane to be in the ground state initially. 



beam splitter interaction. As indicated by this numer- 
ical example, driving the system with time-modulated 
driving strength provides a useful method for generat- 
ing entanglement in a quantum system that is in contact 
with a thermal bath. 

For the parameters discussed here, the rates of the 
optomechanical cooling and heating processes are equal 
(F± ~ Fc), and these processes reduce the entanglement. 
We note, however, that the optomechanical damping can 
in principle also be used to generate entanglement in a 
steady-state situation (i.e. for fixed coupling G) by re- 
ducing the effective occupation numbers of the mechani- 
cal oscillators as discussed in [371 . 



TECHNICAL DETAILS REGARDING THE 
EXPERIMENTAL SETUP 



In this section we discuss technical issues regarding the 
specific setup that we have in mind, namely the optimal 
atom location within the cavity field, and the membrane 
heating due to power absorption. 



A. Optimization of parameters 

Demanding strong atom-membrane coupling, i.e. ful- 
filling the conditions G 3> Fc,Fat,Fni (28), in the end 



boils down to satisfying constraints on the cavity and 
membrane geometry, choosing the proper detunings and 
finding suitable atomic transitions. In the following we 
will therefore go through the set of conditions G ^ 
rc,Fat,F„i and Wat = ^^m (p3| in detail. 



In order to obtain weak cavity-induced decay Fc ^ G, 
we concluded in Sec. |III| that it is necessary to drive the 
cavity far off resonance 



A 2> KjWn 



(31) 
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and to keep at the same time a balanced atom-cavity and 
membrane-cavity coupling g-^ ~ ^at, which is equivalent 
to demanding 

Here we first note that the intracavity amplitude a drops 
out. For simplicity we estimate fi ~ 2r and 9{xa^t) ~ 1. 
In the following we insert the respective definitions of g^ 
(151 and C/q (16 1 and use that lo^ — cki, and find 



kL 



kS 



The difference in zero-point fluctuations between mem- 
brane and atom will give a factor ^m/^at = \prnjM. 
Moreover, on the left side, we introduce the cavity finesse 

^ TTC 

T ^ 

and on the right side the cooperativity parameter, 

This way, the ratio of the coupling constants (jm/^at turns 
into a ratio of the cavity finesse T to the reduced single- 
atom cooperativity (7/<5) C, which must be balanced by 
the mass ratio m/M , 



] V M 




(32) 



The equality ( |32[ ) does not only put a condition for weak 
cavity-induced decay, but will also be useful in the fol- 
lowing to connect the respective parameters related to 
the membrane-cavity and atom-cavity coupling. 

Here we should comment on the dependence of the ra- 
tio T jC on the cavity geometry. At first glance it may 
seem like this ratio can be controlled through the cavity 
length L, with T jC ^ {cj L^i^ jVL^. This is however not 
the case. Keeping in mind that the electric field strength 
is proportional to yJ\/V , with the mode volume V = AL 
and A the beam cross section, the dependence on the cav- 
ity length L in the cooperativity C through the relation 
~ 1/{AL) in fact cancels the length dependence in !F, 
assuming fixed cross section A. We find that the rele- 
vant geometric parameter for this ratio is the beam cross 
section. 



J' 
C 



A. 



Next we require small decoherence due to atomic mo- 
mentum diffusion, Fat/G <C 1, which gives the condition. 



First of all, we use the condition on the coupling con- 
stants ^at — 5m to write the inequality ( 33 1 in terms of 
atomic parameters. 



4 {UoavY 
A 



< V Pel- 



The Lamb-Dicke parameter r] drops out. Furthermore, 
with pe = a'^Uo/5, the intracavity field amplitude a also 
drops out, and what remains is a condition on the coop- 
erativity parameter, 



4k' 



(34) 



which has to be very large, taking into account the con- 
dition (31 1. 



Finally, thermal decoherence depends on the ambient 
temperature T of the membrane. As we will discuss in 



more detail further below (see Sec. VCl, it is reasonable 



to assume that heating of the membrane is in fact caused 
dominantly by absorption of laser power, which depends 
in particular on the thermal link Kth of the membrane to 
its support. The condition of small thermal decoherence 
Fm/G <C 1 reads in detail. 



> 



7n 



27r ojcca^ 



(35) 



First of all, we use the condition on the coupling con- 
stants (7at — 5m to write the inequality ( 35 1 without 
atomic parameters. Using also that Wm = ?i/(2M^^), 
we arrive at, 

(2r)2^^a;2a2 _ M^^7,„ tt hLo.ca' 



L2A 



> 



The amplitude a and the zero-point fiuctuation drop 
out of the inequality. In the same fashion as for the 



condition ( 34 1 , we rewrite the inequality with respect to 
A/k, 



(36) 



with the first factor related to the properties of the mem- 
brane in the cavity, the second factor comparing the ther- 
mal link of the membrane to its natural linewidth, and 
the third factor comparing the energy of a single cavity 
photon to the effective "rest energy" of the membrane. 
Remarkably, this condition is independent of the laser 
power, and the left-hand side depends only on parame- 
ters fixed at fabrication. 

Together, Eqs. ([3l|, ([32]), (|34| and ([36| ensure the set 
of conditions for strong coupling in ( 28 1 . Note that the 
intracavity amplitude a dropped out in all cases. The 
absolute timescale of the system is thus not fixed by 



Eqs. (31 1, (32), (34) and (36), but by the resonance con- 
dition tJin ~ cjat which fixes the cavity amplitude a. 



(33) 



2 2 ^'O 

6 
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The membrane frequency Wm is fixed by construction, 
whereas the intracavity amplitude depends on the laser 
power P, 



2P 



B. Details of ac Stark shift potential 

Let us first discuss briefly how the various require- 
ments on the AC Stark potential generated by the two 
cavity modes can be met. These requirements are as fol- 
lows: (i) Above we have assumed for the atom-cavity 
coupling (7at = Uoar]9{xiit) and for the diffusion rate 

2 

rat — "/^^{xat) that both geometrical factors, 9{xs,t) 

and ^{x,t) = [2- (4/5)u(xat)]/0'(a;at) (for Am = 0), can 
be of order one for a proper choice of the atomic mean 
position Xat along the cavity axis. Moreover, it is desir- 
able to keep the value of which enters the atomic 
trap frequency, as well close to one. (ii) The two modes 
have to couple, respectively, to the Di and D2 lines of 
the chosen atomic species. For a micro-cavity this implies 
that the two modes have to be separated by a couple (say 
q) of free spectral ranges (FSRs) only. A typical inten- 
sity profile is shown in Fig. |9^ for a mode separation of 
q = 5 FSRs. (iii) The atom has to be located in one of 



u{x) 
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FIG. 9; (a) Spatial dependence u{x) of the AC Stark poten- 
tial along the cavity axis for a cavity length L ~ 53/j,m. The 
two driven modes are at Ai ~ 852nm and A2 — 888nm. Their 
separation is q=5 FSRs. Both the parameter 6{xa_t) (b) enter- 
ing the atom-cavity coupling gat and the parameter ^(Sat) (c) 
entering the atomic dephasing rate Fat can be kept close to 
1 at potential wells around the points where Skx — nn with 
Sk = ki — k2 and n < q. In (d) the parameter C^(x) is shown, 
which can be well around 50% while 9, f ~ 1. The resulting 
loss in trap frequency can be easily compensated for by an 
increased intracavity amplitude. 



the potential wells determined by u'{x) = 0. If we take 
^1(2) — {k ± dk)/2 then this condition is equivalent to 



tan(kx) 



Sk 



ta,n(Skx). 



In Fig. [9]d, c and d we show exemplarily the values 
of the parameters 9, ^ and ^ for the possible potential 
wells, i.e. for the solutions to the last equation. As can 
be seen the intensity maxima which exhibit the desired 
properties lie close to points where the cavity modes are 
almost completely out of phase, that is at points where 
Skx = mr for some natural number n < q. 



C. Membrane heating due to laser absorption 

Thermal decoherence depends on the ambient temper- 
ature T of the membrane. It can be reduced by precool- 
ing the membrane with a cryostat. However, it is im- 
portant to note that there is a natural lower limit for T 
which is set by absorption of laser light inside the mem- 
brane. The intracavity light hits the membrane in its 
center, where a fraction a — Pa/Pc of the overall cir- 
culating power Pc = in the two cavity modes is 

absorbed. If the cavity finesse J- is limited by absorption 
inside the membrane, we can estimate a ~ The ab- 



sorbed power Pa causes an increase of the temperature of 
the membrane center by AT ~ ■; — - — P„, where Kt-h is the 
thermal link of the membrane center to the membrane 
supporting frame |39j . Kth depends on the specific geom- 
etry and the material properties and is chosen here such 
as to have dimensions of Hz. While it is not entirely clear 
how the resulting inhomogeneous temperature distribu- 
tion exactly affects the vibrational mode in question, a 
safe assumption is an increase of the ambient tempera- 
ture by AT. 

In |39j , experiments on heat transport inside SiN mem- 
branes at cryogenic temperatures were performed. By 
rescaling the thermal link measured in |39] at a tempera- 
ture of ~ 2 K to our geometry, kB^th — 10 nW/K is ob- 
tained. We furthermore use our parameters Pc = 850 /iW 
and ^ = 2 X 10^, noting that this value for is consis- 
tent with an imaginary part of the refractive index of the 
membrane of Im(rt) ~ 1 x 10~^ [T^]. With these parame- 
ters, we obtain AT ~ 2.5 K. Cryogenic precooling of the 
membrane frame to Tq < AT thus allows one to obtain 
membrane temperatures of the order of T ~ AT. 

To gain further insight into the temperature distribu- 
tion inside the membrane, we simulate heat transport in 
the membrane by solving the heat equation in 2D with 
the finite elements method. We assume that the absorbed 
power Pa is homogeneously distributed over an area 



A = 



in the membrane center, where 



Wo 



10 /im 



is the beam waist of the cavity mode, and that the mem- 
brane frame is held at a fixed temperature of Tq = 2 K. 
At this temperature, the thermal conductivity of SiN is 
kth = 0.05 W/mK [39J. 
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Figure [TO] shows the steady-state temperature distri- 
bution T(y, z) in the membrane obtained from the sim- 
ulation. The peak temperature in the membrane center 
is T(0, 0) = 5.8 K. The average temperature obtained 
by integrating T(y, z) over the membrane cross sectional 
area is T = 2.8 K. We note that our simulation overesti- 
mates the temperature increase, because a constant value 
for fcth is used, while in reality fcth increases rapidly with 
temperature |39j . In summary, we conclude that mem- 
brane heating due to laser absorption sets a lower limit 
on the attainable T, but for our parameters still allows 
for cryogenic precooling to T of a few kelvin. 

Temperature T(y,z) [K] 




FIG. 10: Finite elements simulation of membrane heating due 
to laser absorption. The steady-state temperature distribu- 
tion T{y, z) in the membrane is shown for the experimental 
parameters considered in this paper. The laser spot of radius 
Wq = 10 fim is indicated by the black circle. The membrane 
frame is held at fixed To — 2 K. 



VI. OUTLOOK AND CONCLUSIONS 

In this paper we have discussed the coupling of the 
motion of a single atom and a mesoscopic mechanical os- 
cillator, giving rise to a coupled oscillator dynamics. As 
an outlook, we want to indicate how the present setup 
could in principle also be used for implementing a Jaynes- 
Cummings model by coupling the membrane vibrations 
to the internal atomic degrees of freedom. Consider an 
atom with two stable ground states trapped by an exter- 
nal potential inside the cavity. Let both levels be Stark- 
shifted by the cavity mode, but in opposite directions. 
The coupling of the cavity field quadrature to the two 
level system is then given by 

— a(a -I- a' jaz- 
d 

Changing basis and doing the rotating wave approxi- 
mation, we find an atom-cavity interaction of Jaynes- 



Cummings form, 

9at{a(T+ + aV_), gat = UqU. 

From this brief derivation we learn that by coupling di- 
rectly to the internal levels one wins a Lamb-Dicke pa- 
rameter ?7 <C 1 in the coupling constant, compared to 
coupling to the atomic motion. However, this comes at 
the prize of an increased dissipation rate, as the sponta- 
neous emission from excited states with rate 7 translates 
into ground-state dephasing {Tg^t/2)T>[az]{p) with rate 

Tat ~ IPe 

which is a factor 1/ry^ larger, compared to the momen- 
tum diffusion. In order to profit from the increased cou- 
pling strength, it would thus be crucial to further sup- 
press spontaneous emission by techniques such as e.g. 
coherent population trapping. Our main point here is 
however not to declare a winning model class, but rather 
to point out the various types of interaction which can 
be implemented with the present setup. Realizing the 
Jaynes-Cummings model as described above allows for 
experiments along the lines of microwave CQED [38,. 

Overall, our results illustrate that the present setup 
actually provides a toolbox for engineering various in- 
teractions and different types of dynamics, which pave 
the way towards quantum state engineering of and full 
quantum control over massive micromechanical systems. 
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APPENDIX A: ADIABATIC ELIMINATION OF 
THE CAVITY MODES 

In this section we present details of the adiabatic elim- 
ination of the (independent) cavity modes, starting from 
the linearized master equation in the interaction picture 
w.r.t. Hq, 



W = -i[Hint{t),W] 



and further assuming atom and mirror to be on reso- 
nance, Wat = ^m- The method wc employ here is the 
projection technique |29j which assumes separation of 
timescales; that the fast cavity dynamics on the time 
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scale of the interaction allows us to approximate the den- 



sity operator as W 



p for times t > l/A^. Here 



Pc = 'S>i\0)i{0\i is the steady-state of the shifted cavity 
modes, and p — Trc(Vl^) is the reduced density operator 
for atom and membrane motion. Furthermore, the in- 
fluence of the cavity on the membrane- atom dynamics, 
i.e. the correction in p to the free dynamics, is included 
through second-order perturbation expansion in the in- 
teraction term, using gi/{Ai ± Wm) ^ 1- Finally the 
expressions are simplified using the Born-Markov approx- 
imation. A crucial point is the assumption that the cav- 
ity dynamics dominates over the independent dissipation 
mechanisms of membrane and atom, lAzbtJml 3> FatiFm, 
which would otherwise complicate the cavity-mediated 
dynamics considerably. The result of the adiabatic elim- 
ination is an effective master equation for the atom- 
membrane system, 

P = {Liu + L^t + ic-mod) (p), 

with cavity-mediated atom-membrane dynamics, 
1*00 

Jo 

{ [Hiutit), {®^e'^'^) [Hiutit - r),p° ® p{t)]]} . 

After tracing over the cavity modes, we find that due to 
the independence of the cavity modes (only combinations 
of the form Trc{ajaip[!} contribute) the cavity-mediated 
Liouvillian is a sum of contributions from the different 
modes «, 



+ F^t)) , Pit) [E,{t -r)+ F}{t - t) 



Performing the integrals and returning to the lab frame, 
we find 



1) P {h-,^F^ 



h.c. 



[p{ 



h-_,F, 



h+^F^ 



F, + F 



h.c. 



with A^J2^ {h^,^F^ + h+^,F}^ (f, + f}^ . Reading off 

from the commutator that iJat-^m — ~ ^^); 

find the cavity-mediated coherent dynamics. 



h^,F, + h+,iF,, 



F 



h.c. 

(A3) 



The anti-commutator from the second line of (Al) com 



bined with the sandwich terms in the first line of (All 



describes correlated decay of membrane and atomic mo- 
tion through the cavity. 



2 F, 



,F, + h+^,Fl 



f}^ P (h^,,F, 4 



h+,^F] 



F, + F, 



h.c. 



(A4) 



When written in this form the decay is not manifestly on 
Lindblad form I?[a](/9) = 2apa^ — {a^a,p}+. However, it 
is possible to diagonalize the Liouvillian (A4| and write 



it as a sum of 2 independent jump processes {k = 1,2) 
for each cavity mode i, 



Lc{p)^J2^'^ 



with the jump operators 



(0 



J, 



F, 
F} 



which are described by the eigenvectors m)^, of the 2x2 
matrices M,-, 



A/,; 



2Re{/i+ -I- h^S, 
{h^^, + hX^,)* 2Re{;i_,,} 



(A5) 



The corresponding decay rates 7^ are given by the eigen- 
values of Mi. However in view of the fast rotations of 
terms of the form F^pF^ for ld-ui 3> 7^*\ we perform the 
rotating wave approximation (RWA), which here corre- 
sponds to omitting the off-diagonal terms in Mi, leading 
to (p5|. 



(Al) 



with the constants h± i ~ /[i^ + i(Aj ± 10^^)] ■ 

We expect the cavity-mediated dynamics to be de- 
scribed by 



ic-mod(p) = ~i [i?at-m,p] + (p) , 



(A2) 



with cavity-mediated interaction i/at-m (including cor- 
rections to the free dynamics) and cavity-mediated decay 
Lc{p)- In order to compare the expectation ( |A2[ ) with the 
result ( Al ), we split the second line of ( Al I into commu- 



tator and anticommutator parts. 



APPENDIX B: TIME EVOLUTION AND STATE 
TRANSFER 

The reduced atom-membrane master equation Q can 
be written on the general form 



R'^HR, p 



ip) 



(Bl) 



(pA + h.c.) 



{A-A^),p 



{A + A^),p 



with the matrix H describing the Hamiltonian dynamics 
and the vectors describing the jump operators of the 
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dissipation processes. At this point it is convenient to 
switch to the dimensionless {X, P} language, and thereby 
introduce the basis vector R = [-^m, ^m, -''^at, -Pat]"^ with 
Xat = (flat + aIt)/^/2 and ^ (a,„ + al^)/V2 and the 
commutation relations collected in a matrix a, 





' 1 








-1 




1 


iaij, a = 




1 








-1 0. 


The 



The time evolution of a Gaussian state is fully described 
by its covariance matrix j(t) and displacement vector 
d(t), which are defined as the first and second moments, 
respectively. 



Hence dissipation does not enter the matrix Q ( B2 1 . Con- 
sequently the time evolution of the displacement vector 
d(t) is completely coherent, and the effect of dissipation 
on the state can only be seen in the evolution of the co- 
variance matrix ^{t). 



Transfer of coherent or squeezed states 



atom-membrane interaction is here studied by solving the 



equations of motion (B3) for the covariance matrix ^(t). 



d= (R), 



%■ = ^(R.,Rj 



For an initial coherent state, 



7(0) 




From the general master equation ( |B1| one can derive we find the following evolution 



the equations of motion for the moments, 

3(t) = Qd(t), 4(t)=Q7(t) + 7(t)Q^ 
with the matrices Q and N given by. 



■N, 



2a (H + Im{r} 



r 



N = 2(7 



Re{r} 



(B2) 



Here the matrix F collects the information about the var- 
ious dissipation channels, 

Tfc 
^ 2 

k 



Using this technique, we can solve for the time evolution 
of Gaussian states. In particular the solution for the time 
evolution of the covariance matrix reads. 



We now go on to present the analytical results for 
state transfer which will be used in section |IV| based 
on the calculations in App. [B] Along the lines of sec- 
tion |IV| we consider the symmetric two-mode setup with 
5m = ffat (= g/V^), assuming the rotating wave approxi- 
mation Will ^ G, Fat,7fc*'' and the large detuning regime, 
|A| 3> (WmjK, and neglecting the difference in membrane 
heating and cooling, 1/nm ^ 1. For this special case, 
the jump operators form cooling/heating pairs with equal 
rate F; if a cooling channel is described by some vector 
L, then the corresponding heating channel is described 

by its complex conjugate ^L^ , 



P' 



F 
2 



V 



{p)+V 



R 



Thus, in this particular case, the sum of the respective 
contributions to the dissipation matrix F from the cooling 
and the heating processes is by definition real, 

Im{f } = 0. 



7(i)--[l + 2n(i)] 



+ 7corr(t) 



describing a thermal state with dissipation-induced pop- 
ulation 

n{t) = ^ (2Fe + F,n + Fat) t. 

The second part of the covariance matrix, 7corr , describes 
oscillating atom-membrane correlations which are only 
present when atom and membrane dissipate with differ- 
ent rates. 



7(t) = eQ*7(0)eQ"*+ / dTeQ(*-^)NeQ"(*-^). (B3) % 





( si{t) 




S2{t) \ 


Tm — Tat 




si{t) 




2G 




-S2{t) 


-Slit) 











with si{t) = sinpGi] and S2{t) = sin^[G't]. Due to ther- 
malization the fidelity of the state transfer decreases with 
time from F(0) = 1 according to 



v/det (7mW + 7at(0)) 



(B4) 



as presented for a state swap (t = ■k/2G) in Eq. (30 1 
and Fig. |4]jc). Note that this measure does not take 
into account the (coherent) rotation of the displacement 
vector. 

For the transfer of a squeezed state we see a similar 
thermalization, as is clear from Fig. [5] where we study 
swap of a state with an initial atomic quadrature squeez- 
ing of e^^ 



7at (0) = 
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In this case the covariance matrix evolves according to 

/ C ~t~ 2 /Is wan 1 



2ns 



(B5) 

Here 7,'ji(i) is the diagonahzed membrane covariance ma- 
trix, and the thermal population is given by 



''swap 



^"i(^) - ^(2rc + r^ + r,o. (B6) 



2. Transfer of a Fock state 

Our figure of merit for the thermalization of a Fock 
state during state transfer to the membrane, is the neg- 
ativity of the membrane Wigner function Wm(/3,/3*,t) at 
the origin /3 = /?* = relative to the (absolute) negativ- 
ity of the Fock state Wigner function u;f(0, 0), 



iVw(i) 



hF(0,0)| 



1 

47r 



iU.t). (B7) 



Here Xm(Cm, i) is the characteristic function of the mem- 
brane state, which is derived from the characteristic func- 
tion of the total atom-membrane state, 

Xm(Cm,t) = Xm.at('?m,'?at = 0,i). 

In order to evaluate this expression we need the time evo- 
lution of the full characteristic function, which is defined 
in terms of the density operator /Oni,at(^), 

Xni.at (6n,^at,^) = Tr{p,n,at(<)£'m(Cm)£'at(^at)} 

with -D(^) = exp [^a^ — ^*a] the displacement operator. 
Let us first consider the coherent membrane-atom state, 
for which the time evolution is described by the covari- 
ance matrix ^{t) and the displacement vector d(i). 



Xa,/3('?m,Cat,t) = exp 



with 



( Re(e„,) \ 
Ini(em) 

V I m(£at) / 



Our next attempt is to use the result for coherent states 
to simplify the description of the evolution of a Fock 
state; here we focus on the state |1). We note that by 
expressing the Fock state in terms of coherent states. 



\l)^dje\-\''^\a)) , 



the characteristic function Xi(^) for |1) can be written in 
terms of corresponding functions for coherent states, 

xm - dl^, (e""*Xa(0) ■ (B8) 
Consequently, for an initial density matrix pyS) — 
|0)(0|in|l)(l|at we derive the initial characteristic func- 
tion x(Cm,Cat,0), 



X(Cm,Cat,0) 



d^a* I e"" exp 
with initial conditions 



7(0) = 2 



(B9) 














d(0) = V2 



( \ 



Re(a) 
V Ini(a) ) 



Using the linear properties of the time evolution map et, 

p{t) = 6t(p(0)) 

we find that the time evolution of the full characteristic 
function is given by 



X(fm,^at,t) 



di 



exp 



a=0 



from which we derive the membrane Wigner function 
negativity (B7); for this specific case it reads. 



2n{t) + cospGt] + sin[2Gt](r„i - rat)/2G 
(1 + 2n{t) + sin[2Gt](r^ - rat)/2G)' 
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